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We theoretically demonstrate that atomically-precise 
"nanoscale processing" can be reproducibly performed by 
adaptive laser pulses. We present the new approach on the 
controlled welding of crossed carbon nanotubes, giving various 
rnetastable junctions of interest. Adaptive laser pulses could 
be also used in preparation of other hybrid nanostructures. 
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Understanding and fully controlling growth, structural 
modifications and coalescence of nanoscale materials has 
a top technological priority. Many excellent examples of 
such systems can be found among the numerous types of 
recently discovered nanotubes [1-3] . These unique mate- 
rials are usually grown under rather poorly understood 
conditions using plasma [4], laser [5] and other strongly 
non-equilibrium processes, where catalytic atoms can ac- 
tivate "zipping" of chemical bonds [6]. 

Although artificial nanoscale materials resemble bio- 
macromolecules, they do not possess "inheritable mem- 
ories" , analogous to DNA, crucial for their exact repro- 
duction. Thus our principal question is, if we can develop 
effective atomically-precise processing methods, that can 
reproducibly prepare such systems. This goal has been, 
for example, achieved in a molecular beam growth of su- 
perlattices, where nanoscale patterns in one dimension 
are practically under full control [7]. 

It would be especially attractive to achieve a precise 
control over some advanced structural modifications. An 
example is coalescent "welding" of tubular structures [8] , 
which could lead in crossed nanotubes [9] to prepara- 
tion of strong light-weight nets. As just demonstrated 
experimentally [10], atomically-smooth welded nanotubc 
structures can be induced by irradiation with electrons 
of MeV energies. Unfortunately, control over this process 
is so far limited, and thermal healing of the radiation- 
induced defects [11] is only partial. 

It thus becomes very interesting to consider thermal 
welding of nanotubes, driven by an externally applied 
pressure. Recent simulations [12,13] have revealed that 
this approach could produce defect-free welded struc- 
tures, with deep energy minima. On the other hand, 
practical applications could also largely benefit from a re- 
producible preparation of various rnetastable atomically- 
precise structures, with shallow energy minima. An ex- 
ample are nanotube junctions with a "quasi-continuous" 
variation of their atomic structures. 

Here, we address this challenging task, and explore the 
possibility of a reproducible welding of nanotubes, by 
"adaptive" laser pulses. A feedback control of ultrashort 




FIG. 1. Coalescent welding of crossed single-wall (SWNT) 
or multi-wall (MWNT) carbon nanotubes induced by adap- 
tive laser pulses. Local time-dependent excitation heats a 
micron-size region of the nanotubes, and enables controlled 
atomically-precise restructuralization of C bonds at the pres- 
sured tube crossing (thin dashed lines). 

chirped pulses is already applied in traditional macro- 
welding of thin films [14]. A precise nanoscale processing 
could, in principle, be performed by more sophisticated 
pulses [15], prepared in optimal control techniques [16], 
that can selectively break molecular bonds in gases and 
liquids [17]. Such pulses possess many degrees of freedom 
in their complex shaping and tuning, with the potential 
for a storage of "production informations" . 

In Fig. 1, we present a scheme of the controlled welding 
of crossed nanotubes, realized by their local excitation 
with adaptive laser pulses. These pulses induce interband 
electron transitions in a micron-size region, where the hot 
carriers emit LO phonons. The generated phonons lo- 
cally decay and heat nanotubes below their melting point 
(Tmeit ~ 4000 K) [18]. Reconstruction of C bonds is 
thus induced thermally in the contact area, where nu- 
cleation energies of potential defects are decreased by 
external pressure. Since the flipped bonds are selected 
and directed by the configuration of the whole system 
and the applied pressure, light pulses mostly control the 
total extent of coalescence. The joined region preserve 
the hexagonal atomic pattern if the tube structures are 
complementary (armchair/zigzag). 

We start the analysis of the welding control with the 
description of the nanotube excitation by inhomogeneous 
light intensity E(x,t) [19]. From the Fermi's Golden 
rule we can estimate the electron/hole injection rate, 
h e (h)(x,t) oc \E(x,t)\ , in the conduction/ valence band 
of the crossed SWNT. Typically we use the field inten- 
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sity Eq — 5 MV/cm at frequency u> = 10 15 s -1 , for which 
we obtain h e (h) ~ 0.3 ps -1 per unit cell, with 40 C atoms 
for the (10,10) nanotube [20]. This relatively strong exci- 
tation could cause variation of the optical coupling, and 
even lead to a dielectric breakdown of the system. Never- 
theless, it should be possible to adapt the light pulses to 
these conditions in the learning process, so the structures 
could be controllably prepared without destruction. 

The hot electrons and holes are generated with energies 
of a few (I « 3) LO phonons, which are emitted within 
t « 0.5 ps. Since the carrier velocity is v e « l^m/ps, 
their energy is relaxed within the micron-size generation 
volume. Therefore, the phonon injection rate n p (x,t) 
approximately follows the excited electron/hole profile 
h p (x, t) w I (n e {x, t) + hh(x, t)) . The LO phonons decay 
locally into LA phonons, which carry most of the heat. 
Thus the linear density of heat H(x, t), generated in each 
of the crossed nanotubes, is of the form 

Tt(x, t) « £ p h p (x, t) oc 2l£ p \E(x,t)\ 2 , (1) 

where £ p = Iivq is the LO phonon energy. 

We can describe the one dimensional transport of heat 
by the equation 



l( K{ T)^) + H(x,t) = C { T )p 



dT 



(2) 



Here, k(T) is thermal conductivity, C(T) specific heat, 
S w 1.18 nm 2 the effective surface per tube in a rope and 
p = 1.9 • 10~ 23 kg/nm the mass density of the considered 
(10,10) armchair nanotube. Both theoretical [21] and ex- 
perimental [22] results in nanotubes lead to a large k(T), 
dominated by a phonon transport. From the predicted 
value k(T = 300 K) « 7500 W/K-m in SWNT, adjusted 
to the high-temperature profile in graphite [23], we ob- 
tain k(T) w 23.5 • 10 3 /T 1/5 W/K-m. We can also use the 
specific heat of graphite [23] , with the high-temperature 
fit C(T) w 830 + 28.2 (T - 300) 1/2 J/kg-K. 

The thermally induced coalescence of nanotubes is 
driven by lessening of the pressure-assisted potential en- 
ergy of the system (enthalpy), as in sintering. It can be 
realized by sequences of created and annihilated Stone- 
Wales (SW) 5/7 defects [13], and some other types of 
defects [10]. Typically, the energy of the SW defect is 
Ef w 2 — e Co eV, where e is a local strain [12]. The for- 
mation of defects is limited by the size of their activation 
barrier E* w 6 — e Cq eV, which can block the process 
even if Ef < 0. The last condition can be met already at 
strains e « 0.1, since |C | ~ \C$\ « 20-30 eV. Numerical 
simulations show [12] that the magnitude and sign of the 
material parameters Co and Cq depend, in general, on 
the defect orientation and the tube chirality. Therefore, 
at the tube crossing, where the pressure-induced strain e 
is large, only defects with large and negative Co and Cq 
would be thermally nucleated. 

This fact could significantly reduce the large number 
of possible coalescence paths, switching between various 



intermediate structures. The nucleated defects are pre- 
dominantly chosen by the types of nanotubes, geometry 
of their crossing and the applied pressure. Therefore, 
fluctuations of the followed coalescence paths should be 
relatively small, so that welded structures with the same 
number of defects could be reasonably similar. The most 
favorable path contains structures with the lowest ener- 
gies. For example, a head-to-head coalescence of two 
(10,10) C nanotubes can be realized by following a path 
with just 68 steps of SW defects [13]. In this process, 
the flipped bonds first interconnect the tubes, soon after 
a neck is formed, and finally the tubes become smoothly 
rejoined (see Fig. 1). The goal of the optimal control is to 
deliberately stop this process at any of the metastable in- 
termediate structures, formed in the vicinity of the most 
favorable path. 

Consider that the system occupies with the probability 
Pi any of the structures with i defects, formed in the 
welding process. The probability pi is described by the 
rate equation 



Pi-i , Pi+i 



i+l 
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where the nucleation time r 4 + (rf) gives the transition 
rate for creation of i + 1 (i — 1) defects, if the system has 
i defects. The nucleation processes are activated by the 
temperature at the tube crossing, so the times are 



(rf)" 1 « 2N at v Q exp [- E* ± (i) / k b T] . 



(4) 



Here, N at w 20 is the effective number of C atoms in the 
local region, vq = 10 ps -1 is the vibrational frequency, 
used also in £ p , and E±(i) are activation barriers of the 
i + 1-th and i — 1-th defects for structures with i defects. 
We model these barriers as follows, 



£?i(i) = 6-Aexp[(-|i-io|±0.5)/A] eV, 



(5) 



where we assume that the pressure- induced barrier shift, 
A > 0, is exponentially relaxed with the growing num- 
ber i of defects, since the pressure (motion of the tubes) 
cannot follow the fast (pulsed) coalescence. In numer- 
ical testing, we apply ns welding pulses, separated by 
/is periods. After the first pulse, the strain e can re- 
turn to its initial value in all structures, irrespective of 
their number of created defects jo- ln the next pulse, 
structures with different number of defects jo thus ef- 
fectively start their evolution with the same barriers 
E* ± {i = j ) = 6 - Aexp[±0.5/A]- Here, the shifts ±0.5 
reflect the pressure-induced asymmetry in transitions in- 
creasing/decreasing the number of defects by one, which 
is the driving force of the coalescence. 

Let us now discuss in more details practical realization 
of the welding control. In experiments, we can search the 
optimal field E(t) in a feedback learning loop (see inset 
of Fig. 4), where we impose our welding requirements. 
In contrast to the MBE [7], where growth of individual 
monolayers is controlled in situ, here the light pulses need 
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FIG. 2. Time-dependent temperature T(x,t) at the cross- 
ing (x = 0), induced by an optimal pulse. Three such (sepa- 
rated) pulses create in total (TV) ~ 20 defects. Spatial distri- 
butions of T(x,t) at different times are shown in the inset. 

to be adjusted to the experimentally obtained structures. 
Mathematically, this search minimizes certain functional, 
which reflects the welding requirements [16] and enables 
theoretical testing of the method. Since our goal is to 
reproducibly obtain certain metastable structures, i.e. 
structures with approximately the same given number 
Nd of defects, we look for a field E(x, t) which minimizes 
the difference 

F = \(Af)-Af D \, (^)=J2 Pi i n . (6) 

i 

Here (J\f n ) arc moments of the distribution pi. The vari- 
ance o~d = \J (A/" 2 ) — (TV) 2 describes fluctuations of de- 
fects in the prepared structures. 

In order to control the nanotube welding, we can take 
advantage of the fact that the pressure-induced lowering 
of the barriers E±(i) stops during the fast creation of the 
first Di defects. Thus, we use short light pulses, heating 
the system in such a way, that less than Di defects are 
created by each of them, so that fluctuations in <jd re- 
main suppressed by the growing barriers E±(i). If we try 
to induce more than Di defects by a single pulse, then the 
tubes would have to be largely heated, since the barriers 
are not lowered any more by pressure. Therefore, struc- 
tural fluctuations would significantly grow and welding 
could spread without control in the whole region of the 
crossed tubes. This means that metastable structures 
with more than Di defects would have to be prepared by 
several pulses, each generating (J\f) < Di defects. 

We now test these ideas by solving Eqns. (1-6). Since 
the flipped bonds are largely selected by the system con- 
figuration (see above), the field practically only controls 
the progress of welding. Therefore, we can search it in a 
Gaussian form E(x,t) = E cxp(~x 2 /er 2 - (t-t a ) 2 /af), 
where we fix the parameters cr t — 2.5 ns (to = 4 ns) 




FIG. 3. Time-dependent probabilities Pi(t) of structures 
with i defects, obtained by excitation with three ns pulses, 
such as in Fig. 2. The pulses are separated by ps periods, 
denoted by the dotted lines. With every pulse, the average 
number of defects moves up, until it reaches (TV") ~ 20 defects. 

and o~ x = 1 /im, and thus choose Eq as the only varied 
parameter. In each iteration, of the learning scheme in 
Fig. 4, we apply three identical (separated) pulses to the 
system. Then we calculate the difference F, and let Eq 
to evolve until F is minimized. In real experiments, we 
can also vary the shapes of the pulses, in order to fit more 
complex transient conditions. 

Here, we assume that the nanotubes are in a contact 
with each other at the crossing and with four contacts 
(two for each), separated 3 /xm away from the tube cross- 
ing. They pull the two tubes in opposite directions and 
thus maintain the vertical force F — 5 — 15 nN [13] at 
their crossing. We assume that the force leads to the 
model parameters A w 4 eV and Di « 15, used in 
Eqn. (5). The contacts also provide heat sinks, and are 
thus held at the temperature To = 300 K. In this geo- 
metrical configuration, the nanotubes move and restore 
the force on /is timescale, after each welding pulse. 

Figure 2 shows the temperatures T(x = 0, t) at the 
crossing, induced by excitation with any of the three (sep- 
arated) laser pulses. By requesting that the pulses create 
in total (./V) « 20 defects, we have obtained their com- 
mon optimal field amplitude, Eq rs 5 MV/cm, in several 
tens of iterations. In the present pulsed regime the tem- 
perature rises up to T w 3000 K, then it slowly relaxes, 
as the heat diffuses through the tubes, and is pumped out 
through the contacts. In the inset, we show the related 
broadening and reshaping of the temperature profile. 

In Fig. 3, we show the probabilities Pi(t) of prepar- 
ing structures with i defects, by application of the above 
three pulses separated by /xs time periods. The pressure- 
induced coalescence leads to narrow population maxima 
at certain numbers of defects i after each pulse. The pulse 
creates less than Di defects, before the pressure is tem- 
porarily released and the coalescence stops. When the 
pressure rises again, the process continues with the next 
pulse, until (A/") « 20 defects arc created. In Fig. 4, we 
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FIG. 4. Time-dependence of the average number of defects 
(A/ - ) (solid line) and their variance an (dashed line) for the 
pulsed excitation in Fig. 3. In the inset, we show the learning 
scheme used in optimal control of the welding process. 



also present the time-dependence of the average number 
of defects (AT) and the variance an, for the excitation 
in Fig. 3. The simulation shows that (TV) largely de- 
pends on the transferred energy, while an remains small, 
if (AT) < D\ within one pulse. 

In order to speed up formation of nanotube nets, we 
can consider simultaneous irradiation of many such junc- 
tions under a mask, which would concentrate the light 
on the selected spots. Tubular nets might be also pre- 
pared by self-assembly processes, if we succeed to initiate 
nanotube branching, i.e. splitting of their growth in or- 
thogonal directions. 

With more sophisticated light pulses [15-17] the 
method can also selectively catalyze chemical reactions 
inside or on the surfaces of nanotubes. This can lead 
to an efficient preparation of new hybrid materials, with 
on-site grown "filling" and "dressing" . An example is 
light-induced polymerization or growth of nanocrystals 
inside nanotubes [24]. The adaptive pulses could even 
control electronic processes in nanodevices, especially if 
the chemical constituents can be reversibly switched by 
light [25]. 
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